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Abstract 


Nanostructures  with  multiple  resonances  can  exhibit  a  suppressed  or  even  com¬ 
pletely  eliminated  scattering  of  light,  called  a  scattering  dark  state.  We  describe  this 
phenomenon  with  a  general  treatment  of  light  scattering  from  a  multi-resonant  nanos¬ 
tructure  that  is  spherical  or  non-spherical  but  subwavelength  in  size.  With  multiple 
resonances  in  the  same  channel  ( i.e .  same  angular  momentum  and  polarization),  co¬ 
herent  interference  always  leads  to  scattering  dark  states  in  the  low- absorption  limit, 
regardless  of  the  system  details.  The  coupling  between  resonances  is  inevitable  and  can 
be  interpreted  as  arising  from  far-held  or  near-held.  This  is  a  realization  of  coupled- 
resonator-induced  transparency  in  the  context  of  light  scattering,  which  is  related  to 
but  different  from  Fano  resonances.  Explicit  examples  are  given  to  illustrate  these 
concepts. 

KEYWORDS:  Scattering  dark  states,  nanostructures,  light  scattering,  coupled-resonator-induced 
transparency,  Fano  resonances,  temporal  coupled-mode  theory 
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When  macroscopic  structures  are  shrunk  to  the  nanoscale,  their  optical  properties  depart 
dramatically  from  the  intuitive  ray-optics  picture. 1  Subwavelength  structures  on  resonance 
can  have  scattering  cross  sections  much  larger  than  their  geometrical  sizes,2,3  and  the  pres¬ 
ence  of  multiple  resonances  leads  to  even  more  possibilities  through  mode  hybridization4 
and  interference  effects.5-9  A  particularly  interesting  phenomenon  is  the  suppressed  scat¬ 
tering  in  nanostructures  with  multiple  plasmonic  resonances, 10-23  plasmonic  and  excitonic 
resonances, 24-30  or  dielectric  resonances,31,32  referred  to  collectively  as  a  “scattering  dark 
state.”  A  wealth  of  models  has  been  employed  to  describe  this  suppressed  scattering,  rang¬ 
ing  from  perturbative  models,12  generalization  of  the  Fano  formula, 13-15  and  electrostatic 
approximation, 22,23  to  coupled-mechanical-oscillator  models. 17  21  These  models  reveal  valu¬ 
able  insights  and  facilitate  the  design  of  specific  structures  with  desired  line  shapes.  However, 
the  general  criteria  for  observing  such  scattering  dark  states  remain  unclear.  Non-scattering 
states  have  been  known  in  atomic  physics  since  the  early  works  of  Fano33  and  have  been 
discovered  in  a  variety  of  nanoscale  systems  in  recent  years5-8,34,35.  However,  Fano  reso¬ 
nances  generally  concern  the  interference  between  a  narrow  discrete  resonance  and  a  broad 
resonance  or  continuum.  Meanwhile,  many  occurrences  of  the  scattering  dark  state  involve 
the  interference  between  multiple  narrow  discrete  resonances,  and  it  seems  necessary  to 
treat  the  multiple  resonances  at  equal  footing.  Thus,  we  seek  a  formalism  analogous  to  the 
phenomenon  of  coupled-resonator- induced  transparency5  that  has  been  established  in  cer¬ 
tain  other  systems  such  as  coupled  mechanical  oscillators36,37,  coupled  cavities38,39,  coupled 
microring  resonators40-44,  and  planar  metamaterials45-48. 

Here,  we  derive  the  general  equations  governing  the  resonant  light  scattering  from  a 
spherical  or  a  non-spherical  but  subwavelength  obstacle,  accounting  for  multiple  resonances 
with  low  loss.  Due  to  the  spherical  symmetry  (or  the  small  size)  of  the  obstacle,  different 
channels  of  the  multipole  fields  are  decoupled.  We  find  that  within  each  channel,  n  resonances 
always  lead  to  n  —  1  scattering  dark  states  in  the  low-absorption  limit.  This  universal  result 
is  independent  of  the  radiative  decay  rates  of  the  resonances,  method  of  coupling  (can  be 
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near-field  or  far-field),  nature  of  the  resonances  (can  be  plasmon,  exciton,  whispering-gallery, 
etc .),  number  of  resonances,  which  of  the  multipole,  TE  or  TM  polarization,  and  other  system 
details.  With  different  choices  of  basis,  one  can  interpret  the  scattering  dark  state  as  arising 
from  the  far-held  coupling  of  multiple  radiating  resonances,  or  arising  from  the  near-held 
coupling  of  a  radiating  and  several  non-radiating  resonances.  We  provide  explicit  examples 
using  plasmonic  resonances  and  whispering-gallery  resonances,  showing  zero  scattering  for 
lossless  materials  and  significantly  suppressed  scattering  for  realistic  materials  with  loss. 
We  also  discuss  potential  applications  for  slow  light,  transparent  projection  screen,  and 
wavelength-selective  transmission. 

Consider  a  linearly  polarized  planewave  incident  on  an  obstacle  that  is  spherical  or  non- 
spherical  but  much  smaller  than  the  wavelength  of  the  incident  light.  This  obstacle  can 
have  arbitrary  number  of  layers  and  material  composition;  in  the  small-obstacle  case,  it 
can  also  be  a  cluster  of  particles.  We  start  with  the  general  formalism  for  such  a  scattering 
problem.  Outside  the  obstacle,  the  electric  held  can  be  written  as  E  =  V  x  (r^TE )-{i/k)V  x 
V  x  (d/)tm),  where  r  is  position  from  the  particle  center,  k  is  the  wave  number,  and  is  a 
scalar  function  satisfying  the  Helmholtz  equation. 1  We  use  cr  to  denote  the  two  polarizations: 
transverse-electric  (TE)  or  transverse-magnetic  (TM),  where  the  electric  or  magnetic  held 
is  perpendicular  to  r.  Choose  our  coordinate  such  that  the  x  axis  and  the  z  axis  is  aligned 
with  the  polarization  and  the  propagation  direction  of  the  incident  light  respectively.  The 
incident  planewave  consists  of  all  multipole  terms  with  l  >  0  and  rn  —  1  in  both  TE  and 
TM1,  so  the  general  solution  can  be  written  as 


VvM,</>)  =  fM)  [si,*h(i1)(kr) 

i=i 


+  4<rh<i2)(kr ) 


Pi  {cos  0), 


(1) 


where  /te(0)  =  sin  <f>,  /tm(0)  =  cos  (j),  (or  lif'1)  is  the  spherical  Henkel  function  of  the 
hrst  (or  second)  kind  corresponding  to  outgoing  (or  incoming)  spherical  wave,  and  P/  is  the 
associated  Legendre  polynomial  with  rn  —  1.  The  amplitudes  of  outgoing  and  incoming  waves 
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sfa  are  coefficients  of  the  general  solution.  Here,  each  angular  momentum  and  polarization 
pair  (/,  a)  corresponds  to  a  distinct  “channel.”  The  spherical  symmetry  (or  the  small  size) 
of  the  obstacle  means  that  different  channels  are  decoupled,  so  optical  response  properties 
are  given  by  the  reflection  coefficients  RiiC  =  sja/ s^a.  Energy  conservation  requires  that 
\Ri,a\  <  1  hr  each  channel.  The  total  scattering  cross  section  of  this  particle  is  given  by 

\  2  00 

*»=  8rEE(2i  +  1)|1-^|2>  (2) 

<y  1=1 

where  Am  is  the  wavelength  in  the  surrounding  medium.  Note  that  there  is  no  inter-channel 
interference  here,  unlike  the  intensity  of  the  scattered  light  at  a  specific  angle  (such  as  back 
scattering)  where  different  channels  of  spherical  waves  can  interfere6,49. 

In  each  channel,  the  scattering  from  a  single  low-loss  resonance  is  given  by  a  Lorentzian 
function2,3,34,35.  Therefore,  when  we  consider  multiple  resonances  each  in  a  different  channel, 
the  total  response  will  be  a  sum  of  Lorentzians  with  no  scattering  dark  state;  this  is  illus¬ 
trated  schematically  in  Figure  la.  For  a  scattering  dark  state,  we  need  multiple  overlapping 
resonances  in  the  same  channel,  as  illustrated  in  Figures  lb  and  lc. 

We  start  with  the  simplest  case  of  two  resonances  in  channel  (/,  a),  as  illustrated  in  Figure 
2a.  The  two  resonances  may  be  of  any  nature  (e.g.,  plasmon,  exciton,  whispering-gallery). 
When  the  resonances  have  low  loss,  R^a  has  two  poles  on  the  complex-frequency  plane,  near 
which  the  system  follows  a  simple  set  of  equations  described  by  the  temporal  coupled-mode 
theory  (TCMT). 2,34,35,39,50  As  a  starting  point,  we  “turn  off”  the  radiation  loss  and  absorption 
loss  of  the  resonances,  and  the  resonance  amplitudes  Aj  (j  =  1,2)  evolve  as 


where  cu'  are  the  resonant  frequencies,  and  u/12,  oj'21  are  the  near-held  coupling  strengths. 
We  proceed  by  changing  to  the  basis  Aj  (without  the  primes)  that  diagonalizes  the  matrix 
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(so  U12  is  now  zero),  which  is  always  possible  because  the  matrix  is  Hermitian  by  energy 
conservation.  In  the  new  basis,  we  “turn  on”  the  low  losses  of  the  resonances,  and  couple 
them  to  the  spherical  wave  in  channel  ( l,  a ), 
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where  7 j  are  radiative  decay  rates,  £j  are  absorptive  decay  rates,  and  Kj  are  coupling  coef¬ 
ficients  to  the  incoming  wave.  Because  both  resonances  radiate  into  the  same  channel  sfal 
the  radiative  coupling  rate  712  is  necessary  to  ensure  energy  conservation.39  Absorptions  do 
not  have  such  constraints,  so  we  let  £12  =  0  for  simplicity.  Meanwhile,  the  outgoing  wave  is 
given  by 

Sl,a  =  4”  d\Ai  +  d2A2:  (5) 


where  dj  are  coupling  coefficients,  and  the  first  term  s~[a  =  sz+CT  comes  from  the  planewave 
itself. 1  Since  we  are  considering  nanostructures,  we  exclude  the  direct  (non-resonant)  back¬ 
ground  scattering  process  related  to  Fano  resonances34,35;  only  the  resonant  scattering  pro¬ 
cesses  are  considered  here.  The  normalizations  are  chosen  such  that  \Aj\2  is  the  electro¬ 
magnetic  energy  in  each  resonance  and  sfa  | 2  is  the  incoming  and  outgoing  power.  For  low 
absorption  loss,  we  can  apply  energy  conservation  and  time  reversal  symmetry39  to  eliminate 
many  unknowns  and  get  Kj  =  dj:  7 j  =  —  d2/2,  and  712  =  —did2/2.  Solving  Eqs  4  and  5  for 
a  steady-state  solution  at  frequency  u  gives  the  reflection  coefficient  s~[a/ s^a  as 


Rlc  =  1 


U(t.  u 


[i(uj2  ~uj)  +  6]  7i  +  [2(^1  -  co)  +  £1] 

\,  i  r  ■  /  \  ,  , 


72 


(6) 


which  yields  the  scattering  cross  sections  spectrum  through  Eq  2.  When  there  is  negligible 
absorption  loss  (£1  =  £2  —  0),  the  particle  becomes  transparent  in  this  channel  (R^a  =  1)  at 
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the  “transparency  frequency,”  defined  as 


CU172  +  CU271 

u ;  =  - 

7i  +  72 

Since  ut  is  simply  a  weighted  average  of  the  two  resonant  frequencies,  it  always  exists, 
regardless  of  the  radiative  decay  rates  and  other  system  details.  We  therefore  conclude  that 
the  scattering  dark  state  is  a  general  phenomenon  in  the  low-absorption  limit  when  two 
resonances  in  the  same  channel  are  simultaneously  excited. 

In  the  presence  of  material  loss,  the  scattering  cross  section  cannot  be  strictly  zero  even 
at  ojf  This  can  be  viewed  as  a  consequence  of  the  optical  theorem1,  which  relates  the  forward 
scattering  amplitude  to  the  total  extinction  cross  section  (scattering  plus  absorption).  Taking 
£1  =  £2  =  £  and  71  =  72  =  7,  the  reflection  coefficient  at  uit  is  R^a  ~  1  —  167£/(o;i  —  u^)2  to 
leading  order  of  the  absorptive  decay  rate  £.  So,  the  “low-absorption  limit”  can  be  quantified 
as  £  <C  (cui  -  U2)2/^. 

The  scattering  dark  state  occurs  when  the  two  outgoing  waves  d\A\  and  ^2^2  have  equal 
magnitude  and  opposite  phase,  as  illustrated  in  Figure  2b.  When  u  <  uq  <  uq  or  when 
uj  >  UJ2  >  uji,  the  two  outgoing  waves  add  up  in  phase  (shown  in  blue  dotted  arrows). 
When  uq  <  u  <  u>2,  the  two  add  up  out  of  phase  (shown  in  green  dashed  arrows).  At 
the  transparency  frequency  cat,  the  two  waves  exactly  cancel  each  other  (shown  in  red  solid 
arrows),  making  the  particle  completely  transparent  in  the  (l,  a)  channel. 

We  now  address  the  choice  of  basis.  In  the  derivation  above,  we  choose  a  basis  to 
express  the  two  resonances  such  that  uj  12  =  0  (no  coupling  in  the  absence  of  loss),  and  the 
scattering  dark  state  arises  from  interference  in  the  far-held  radiation  of  the  two  resonances. 
The  underlying  system  is  independent  of  the  basis,  and  we  may  as  well  choose  a  basis  that 
diagonalizes  the  radiative-decay-rate  matrix  (so  that  7(2  =  0,  no  radiative  coupling;  the 
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prime  denotes  variables  in  that  basis);  in  such  a  basis,  Eqs  4  and  5  become 
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Sl,o  -  St<j  + 

with  cui  =  (cci7i+a;272)/(7i+72),  ^2  =  (^i72+^27i)/(7i+72),  w'12  =  (^-^1)^7172/(71+72), 
7J  =  71+72,  and  k\  =  d\  =  \J d\  +  d%.  The  transformation  of  the  absorptive  decay  rates  is 
similar.  We  see  that  72  =  0  (and  so  k'2  =  d'2  =  0).  So,  in  this  basis,  resonance  A'x  radiates, 
but  resonance  A'2  does  not;  this  is  exactly  the  subradiant-superradiant  model 17,19,21’47’48. 
Intuitively,  this  choice  of  basis  is  possible  because  there  is  only  one  channel  of  radiation,  which 
can  be  incorporated  into  just  one  degree  of  freedom.  The  scattering  dark  state  still  exists  in 
the  low-absorption  limit  at  the  frequency  given  by  Eq  7.  But,  now  it  is  the  direct  near-held 
coupling  >jj\2  that  leads  to  transparency.  We  therefore  conclude  that  the  scattering  dark 
state  can  be  interpreted  as  arising  from  the  far-held  coupling  of  two  radiating  resonances,  or 
arising  from  the  near-held  coupling  of  one  radiating  resonance  and  a  non-radiating  resonance. 
Both  interpretations  are  valid;  it  is  only  a  matter  of  basis  choice. 

The  extension  to  more  than  two  resonances  is  described  in  Supporting  Information.  With 
n  resonances  in  the  same  channel  and  with  negligible  absorption,  scattering  in  this  channel 
vanishes  at  frequencies  given  by 


£  —  =  0. 

Uln  ~  UJ 

]  =  1  J 


(9) 


where  ujj  and  7 are  the  resonant  frequencies  and  radiative  decay  rates.  The  left-hand 
side  swings  from  —  00  to  +00  in  each  interval  between  successive  resonant  frequencies.  So, 
there  is  one  transparency  frequency  in  each  of  these  n  —  1  intervals.  For  example,  there  is 
one  transparency  frequency  with  two  resonances,  two  transparency  frequencies  with  three 
resonances,  etc.  Similar  to  the  two-resonance  case,  it  is  possible  to  choose  a  basis  that 
diagonalizes  the  radiative-decay-rate  matrix;  in  such  a  basis,  one  resonance  radiates  while 


the  other  n  —  1  resonances  do  not — again,  because  there  is  only  one  channel  of  radiation. 

This  brings  our  major  conclusion  that  the  scattering  dark  state  arises  generally  with 
multiple  resonances  in  the  same  channel.  It  does  not  matter  how  many  resonances  there  are, 
or  how  narrow  or  broad  each  resonance  is.  Therefore  the  suppressed  scattering  is  more  general 
than  the  discrete-coupled-to-continuum  or  narrow-coupled-to-broad  scenarios  in  standard 
Fano  resonances33,34. 

The  scattering  dark  state  is  similar  to  “bound  state  in  the  continuum”51-53  because  both 
phenomena  arise  from  destructive  interference  of  outgoing  waves.  However,  a  scattering  dark 
state  is  not  a  bound  state.  Eq  4  indicates  that  the  resonance  amplitudes  decay  to  zero  in 
the  absence  of  the  incoming  wave — they  cannot  sustain  oscillation  on  their  own,  even  at  the 
transparency  frequency.  Also,  the  particle  becomes  transparent  at  steady  state,  but  at  the 
transient  stage  when  d\Ai  +  ^2^2  =  0  has  not  been  established,  the  particle  is  opaque. 

Lastly,  we  provide  a  few  explicit  examples  to  illustrate  the  concepts  discussed  above. 
First,  we  consider  a  multi-layer  sphere  (schematically  shown  in  the  inset  of  Figure  3a)  that 
consists  of  concentric  metallic  core,  dielectric  spacer,  and  metallic  shell.  For  the  first  example, 
we  describe  the  metal  layers  using  the  Drude  model  with  negligible  damping  e(uj)  =  1—u^/uj2, 
let  the  dielectric  layer  be  e  =  2.04,  and  let  the  surrounding  medium  be  air.  For  a  particle  of 
size  [ri ,  rq,  ry]  =  [0.005,  0.073,  0.132]AP  (where  Ap  =  27r c/ay,  is  the  plasma  wavelength,  and 
c  is  the  speed  of  light  in  vacuum),  the  exact  scattering  cross  section  is  plotted  as  the  black 
line  in  Figure  3a;  these  data  are  calculated  using  the  Mie  solution  with  the  transfer  matrix 
method.54  The  electric  dipole  channel  TMi  (l  —  1,  a  —  TM)  dominates  in  the  frequency 
range  plotted;  contributions  from  other  channels  are  six  orders  of  magnitude  smaller.  There 
is  no  absorption  loss  in  this  example,  and  as  expected,  the  scattering  cross  section  goes  to 
zero  at  a  transparency  frequency  ujt  in  between  the  two  resonance  peaks.  We  plot  the  steady- 
state  electric-held  profile  at  one  of  the  resonant  frequencies  and  at  ay  in  Figure  3b  and  3c; 
animations  of  these  held  prohles  are  given  in  Supplementary  Movies  SI  and  S2.  The  held 
prohle  confirms  that  the  particle  becomes  invisible  at  ay.  We  note  that  the  scattering  dark 
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state  is  robust:  perturbations  of  parameters  (such  as  the  layer  thicknesses  and  the  refractive 
indices)  only  shift  the  frequency  where  it  occurs,  consistent  with  our  discussion  above. 

Figure  3a  also  shows  the  prediction  from  the  temporal  coupled-mode  theory  (red  dashed 
line).  We  obtain  the  parameters  in  TCMT  without  doing  curve  fitting.  Instead,  we  locate 
the  two  poles  of  the  exact  reflection  coefficient  R  i  (from  the  Mie  solution)  on  the  complex- 
frequency  plane.  The  pole  locations,  0.4411  —  0.002821  and  0.4456  —  0.002971  (in  units  of 
tCp),  yield  the  TCMT  parameters  uq  =  0.4397,  71  =  0.00285,  0J2  =  0.4470,  and  y2  =  0.00294 
when  compared  to  the  denominator  of  Eq  6.  The  TCMT  prediction  agrees  excellently  with 
the  exact  Mie  solution  even  without  fitting. 

A  collection  of  these  nanospheres  with  subwavelength  spacing  can  act  as  an  effective 
medium.  In  general,  the  scattering  properties  of  such  closely  spaced  nanospheres  are  sub¬ 
stantially  different  from  that  of  an  individual  nanosphere  in  free  space.  However,  when  an 
individual  particle  is  transparent,  we  can  use  the  superposition  principle  to  conclude  that 
a  collection  of  particles  is  also  transparent  even  when  closed  spaced.  So,  near  cct,  we  may 
infer  the  properties  of  the  effective  medium  from  an  individual  particle  in  free  space,  using 
the  standard  mixing  formula55.  Near  ujt,  the  real  part  of  the  individual  particle’s  electric 
polarizability  changes  rapidly  (Figure  3d),  so  the  refractive  index  nes  of  the  effective  medium 
changes  rapidly.  This  leads  to  a  suppressed  group  velocity,  vg  =  du  /  dk  =  c  / (nes-\-udnes / du) , 
for  a  wave  packet  propagating  through  this  medium,  analogous  to  the  atomic  version  of  elec- 
tromagnetically  induced  transparency56,57  and  similar  to  the  metamaterial  realization  using 
subwavelength  optical  antennas46-48.  For  example,  for  a  random  collection  of  the  above- 
mentioned  particles  with  concentration  N  =  (At/3)-3  (where  At  =  2nc/ut  is  the  transparency 
wavelength),  the  group  velocity  is  vg  «  c/200  at  transparency. 

Next,  we  provide  examples  using  realistic  materials  with  loss.  Consider  the  same  multi¬ 
layer  sphere  schematically  shown  in  the  inset  of  Figure  3a,  but  with  the  metal  layers  being 
silver  (complex  permittivity  from  experimental  data58)  instead.  Similar  metal-dielectric- 
metal  nanospheres  have  been  synthesized17,59  and  studied  numerically54,60-63  in  prior  works. 
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To  mimic  potential  experimental  condition,  we  consider  such  nanospheres  suspended  in  an 
aqueous  solution  (e  =  1.77).  Figure  4a  shows  the  exact  scattering  cross  section  for  a  particle 
of  size  [rq,  72,  rs]  =  [20,  31,  52]  nm.  This  particle  is  large  enough  that  the  electrostatic  approx¬ 
imation  is  not  appropriate,  but  small  enough  that  the  TM (  channel  dominates.  Since  the 
absorption  and  radiation  loss  is  relatively  large,  the  prediction  from  TCMT  is  not  as  accurate 
in  this  example;  nonetheless,  we  still  observe  a  clear  dip  in  the  scattering  cross  section  that  is 
suppressed  by  more  than  a  factor  of  ten  compared  to  the  on-resonance  values.  Supplementary 
Movies  S3  and  S4  show  animations  of  the  held  propagation  on-resonance  and  at  the  almost- 
transparent  wavelength,  where  the  particle  becomes  almost  invisible.  Such  nanoparticles  can 
be  useful  for  the  application  of  transparent  projection  screens  using  resonant  scattering  of 
nanoparticles.64  In  the  presence  of  absorption  loss,  the  amount  of  scattering  at  tot  does  vary 
with  perturbations.  As  discussed  earlier,  Ricr  ~  1  —  —  a^)2  depends  on  the  absorp¬ 

tive  loss,  radiative  loss,  and  the  frequency  difference  between  the  two  resonances.  However, 
the  dip  still  exists  with  reasonable  perturbations:  for  the  current  structure,  the  scattering 
minimum  remains  at  least  10  times  smaller  than  the  scattering  maximum  with  a  10%  change 
in  any  thickness  or  refractive  index  parameter. 

Another  possible  application  of  the  scattering  dark  state  is  to  selectively  allow  only  a 
narrow  bandwidth  of  light  to  pass  through  a  medium.  In  a  medium  with  dilute  suspension 
of  particles,  the  transmission  is  given  by  the  Beer-Lambert  law  T  =  exp(—crextNl),  where 
<7 ext  =  Csca  +  ^abs  is  the  extinction  cross  section  of  an  individual  particle,  N  is  its  number 
density,  and  /  is  the  path  length  of  the  medium.  Typically,  resonant  nanoparticles  suppress 
transmission  at  a  target  wavelength.  But  with  the  scattering  dark  state,  we  may  enhance 
transmission  at  a  target  wavelength  instead.  Figure  4b  shows  one  example  using  the  same 
type  of  silver-silica-silver  nanosphere,  designed  to  allow  light  with  wavelength  around  600 
nm  to  pass  through. 

Our  theory  derivation  suggests  that  the  scattering  dark  state  is  not  limited  to  any  par¬ 
ticular  type  of  resonance;  nor  is  it  limited  to  any  particular  radiation  channel.  To  illustrate 
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these  points,  we  consider  a  purely  dielectric  nanostructure  that  supports  whispering-gallery 
resonances  in  multiple  angular  momentum  channels.  Specifically,  we  consider  a  nanosphere 
consisting  of  four  concentric  layers  that  alternate  between  a  high-index  dielectric  (e  =  12, 
can  be  silicon  or  gallium  arsenide)  and  a  low-index  dielectric  (e  =  1,  can  be  a  transparent 
aerogel65).  For  a  particle  of  size  [r'i ,  rq,  73,  74]  =  [40,90,150,160]  nm,  Figure  5  shows  the 
exact  scattering  cross  section  divided  into  individual  channels  as  in  Eq  2.  Scattering  dark 
states  can  be  seen  in  the  TEi,  TE2,  and  TE3  channels.  Note  that  due  to  the  spectral  overlap 
of  the  different  channels,  the  particle  may  become  transparent  to  one  particular  channel  but 
not  the  rest;  in  principle,  it  is  possible  to  probe  the  individual  channels  by  preparing  special 
excitation  waves  (for  example,  a  tightly-focused  beam  can  isolate  the  dipole  channels66). 

In  summary,  we  have  presented  an  analytical  treatment  for  the  general  problem  of  light 
scattering  from  a  multi-resonant  nanostructure  that  is  spherical  or  non-spherical  but  sub¬ 
wavelength  in  size.  With  n  resonances  in  the  same  channel,  n  —  1  scattering  dark  states  arise 
when  the  absorption  loss  is  negligible;  this  result  is  independent  of  the  radiative  decay  rates 
and  many  other  system  details.  The  scattering  dark  state  can  be  interpreted  as  arising  from 
far-held  or  near-held  coupling,  depending  on  the  choice  of  basis.  This  theoretical  treatment 
should  improve  the  understanding  of  scattering  dark  states  and  provide  guidance  for  future 
works  on  this  topic. 

It  will  be  interesting  to  consider  perturbations  that  break  the  spherical  symmetry  of  the 
scatterer,  allowing  each  resonance  to  couple  to  multiple  channels.  The  interference  between 
multiple  resonances  in  multiple  channels  may  lead  to  even  richer  phenomena 11,17,39,67. 
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Figure  1:  Schematic  plots  for  light  scattering  when  multiple  resonances  are  excited  simul¬ 
taneously.  (a)  Two  resonances  in  different  radiation  channels  (e.g.,  TMi  and  TEi).  In  this 
case,  radiation  from  the  two  resonances  add  up  incoherently,  and  the  cross  section  is  a  sum  of 
two  Lorentzian  functions.  The  particle  is  opaque  in  between  the  two  resonances  (illustrated 
in  the  inset),  (b)  Two  resonances  in  the  same  radiation  channel.  In  this  case,  radiation  from 
the  two  resonances  add  up  coherently,  giving  rise  to  a  different  line  shape  and  a  scattering 
dark  state  where  the  particle  becomes  transparent  (illustrated  in  the  inset),  (c)  Three  reso¬ 
nances  in  the  same  radiation  channel.  Here,  two  scattering  dark  states  arise  (illustrated  in 
the  inset). 
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Figure  2:  Illustration  of  a  scattering  dark  state  in  a  doubly  resonant  nanostructure,  (a) 
Abstract  diagram  showing  a  nanostructure  where  two  resonances,  A\  and  A2,  in  the  same 
angular  momentum  l  and  polarization  a  channel  are  simultaneously  excited  by  incoming 
wave  s^a,  and  both  radiate  into  outgoing  wave  sfa.  (b)  Schematic  illustrations  showing  the 
complex  reflection  coefficient  =  s^a/ sfa  with  contributions  from  the  two  resonances 
(each  indicated  by  an  arrow).  When  absorption  loss  is  negligible,  always  lands  on  the 
unit  circle,  and  the  two  outgoing  waves  add  up  either  in  phase  or  out  of  phase.  Different 
excitation  frequencies  are  shown  by  arrows  of  different  colors;  at  frequency  04,  the  two 
outgoing  waves  cancel,  giving  rise  to  a  scattering  dark  state  where  the  particle  becomes 
transparent. 
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Figure  3:  Scattering  dark  state  in  a  doubly  resonant  nanosphere  without  absorption  loss,  (a) 
The  particle’s  scattering  cross  section  in  air,  calculated  by  the  Mie  solution  (black  solid  line) 
and  the  temporal  coupled-mode  theory  (red  dashed  line).  The  inset  schematically  shows 
the  particle’s  composition:  concentric  metallic  core,  silica  spacer  (e  =  2.04),  and  metallic 
shell,  with  radii  [ri,r2,r3]  =  [0.005,  0.073,  0.132]  Ap.  The  metal  is  described  by  the  Drude 
model  with  plasma  frequency  (up  and  negligible  damping.  In  the  plotted  range  of  frequency, 
the  cross  section  is  dominated  by  the  TMj  channel  (electric  dipole),  (b,  c)  Steady-state 
electric-held  pattern,  Ke(Ex),  with  an  incident  wave  Einc  =  elkzx  that  is  (b)  at  the  resonant 
frequency  uq  and  (c)  at  the  transparency  frequency  04.  Animations  of  the  held  propagation 
are  shown  in  Supplementary  Movies  SI  and  S2.  (d)  Electric  polarizability  of  this  particle 
in  air.  At  the  transparency  frequency,  the  large  slope  of  Re(ct)  can  give  rise  to  slow  group 
velocity  for  a  wave  packet  propagating  through  a  collection  of  such  particles. 
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Figure  4:  Scattering  dark  state  in  a  doubly  resonant  nanosphere  with  absorption  loss,  (a) 
The  scattering  cross  section  of  a  particle  consisting  of  concentric  silver  core,  silica  spacer 
(e  =  2.04),  and  silver  shell,  with  radii  of  the  layers  being  [ri,r2,ra]  =  [20,31,52]  nm.  The 
cross  section  is  calculated  with  the  Mie  solution  with  the  particle  in  water  (e  =  1.77)  and  the 
complex  permittivity  of  silver  taken  from  experimental  data. 58  Solid  line  shows  the  sum  from 
all  channels,  and  dashed  line  shows  the  contribution  from  the  TMX  channel.  Animations 
of  the  held  propagation  at  607  nm  and  544  nm  are  shown  in  Supplementary  Movies  S3 
and  S4.  (b)  Transmission  spectrum  for  a  dilute  aqueous  solution  of  nanoparticles  of  size 
[ri,r2,rs]  =  [33,55,77]  nm,  with  the  number  density  times  path  length  being  7  x  109  cm-2. 
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Figure  5:  Scattering  dark  states  in  individual  channels  of  a  nanosphere  made  of  dielectrics 
only.  The  particle  consists  of  four  concentric  layers,  with  relative  permittivities  being 
[ei,  e2,  e3,  e4]  =  [1,12,1,12]  and  radii  being  fry,  r2,  r3,  r4]  =  [40,90,150,160]  nm.  (a)-(c) 
The  particle’s  scattering  cross  section  in  air,  separated  into  different  angular  momentum 
and  polarization  channels.  Scattering  dark  states  can  be  seen  in  the  TEi,  TE2,  and  TE3 
channels. 
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Arbitrarily  many  resonances  in  the  same  channel 

In  the  main  text,  we  consider  the  case  of  two  spectrally  overlapping  resonances  in  the  same  channel.  Here,  we 
consider  the  more  general  case  with  arbitrarily  many  overlapping  resonances  in  channel  ( l, a ).  The  temporal  coupled¬ 
mode  theory  equations  is  now  written  in  matrix  form: 


—  =  (-in-r-E)A  +  KTa+a, 

Sl,a  =  St,a  +  DA‘ 


(S.l) 


Similar  to  the  two-resonance  case,  we  choose  the  basis  such  that  fi  is  diagonal  (so  that  the  resonances  are  orthogonal 
when  we  take  the  loss  to  zero),  and  we  ignore  any  direct  (non-resonant)  background  scattering  process.  With  n 
resonances,  the  matrices  are 
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The  meaning  of  these  variables  are  the  same  as  Eqs.  (4)  and  (5)  in  the  main  text.  With  incoming  light  at  frequency 
co,  Eq.  (S.l)  gives  the  steady-state  reflection  coefficient  sffCT/s'la  as 


Ri,*  =  1  +  D  [i(fi  -  coin)  +  T  +  S]"1  Kt, 


(S.3) 


where  /„  is  the  n-by-n  identity  matrix.  The  particle  becomes  transparent  in  this  channel  when  Ri  a  =  1.  Optical 
theorem  suggests  that  zero  scattering  can  only  occur  with  zero  absorption;  therefore,  we  consider  the  small-absorption 
limit,  5  =  0.  Energy  conservation  and  time  reversal  symmetry  require  the  matrices  to  satisfy  K  =  D,  D  =  2T, 
and  D*  =  —D,  as  shown  in  ref  39.  Therefore,  the  transparency  condition  can  be  simplified  to 


D 


z(n  -  coin)  -  -DTD 


-1 


Dt  =  0. 


(S.4) 


By  rewriting  the  matrix  inverse  using  the  Woodbury  matrix  identity,  we  see  that  this  condition  is  satisfied  when 
D(fl  —  coIn)~1DT  =  0,  which  can  be  written  as 


n 


(S.5) 


because  d2  =  — 2q j.  This  is  Eq.  (11)  in  the  main  text.  This  simple  expression  yields  all  the  transparency  frequencies. 
We  can  immediately  observe  that,  as  we  gradually  increase  co,  the  left-hand  side  of  Eq.  (S.5)  swings  from  — oo  to  +oo 
in  each  interval  between  successive  resonant  frequencies.  So,  there  is  one  transparency  frequency  in  each  of  these  n  —  1 
intervals.  For  example,  there  is  one  transparency  frequency  with  two  resonances,  two  transparency  frequencies  with 
three  resonances,  etc. 

Similar  to  the  two-resonance  case,  it  is  possible  to  choose  a  basis  that  diagonalizes  the  radiative-decay-rate  matrix 
T.  The  matrix  T  =  —^DTD  is  a  rank-one  matrix,  so  only  one  of  its  n  eigenvalues  is  non-zero.  This  non-zero  eigenvalue 
^  -ie;=1^2  =  e;=  7j.  Therefore,  in  the  basis  where  T  is  diagonal,  one  resonance  would  radiate  while  the  other 
n—1  resonances  do  not.  Intuitively,  this  is  because  there  is  only  one  channel  of  radiation. 


SI 


